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A NUMERICAL  DYNAMIC  FRACTURE  ANALYSIS 
OF  THREE  WEDGE-LOADED  DCB  SPECIMENS 


by 

A.  S.  Kobayashi*,  S.  Mall**,  Y.  Urabe***  and  A.  F.  Emery* 


SUMMARY 

A dynamic  finite  element  code  is  used  to  compute  the  dynamic  fracture  tough- 
ness and  crack  arrest  stress  intensity  factor  from  experimentally  determined  crack 
velocities  In  three  fracturing  wedge-loaded  double  cantilever  beam  (DCB)  specimens. 
One  experiment  involving  an  Aradite-B  DCB  specimen  by  Kalthoff,  et  al . , and  two 
experiments  involving  Horaal i te-100  DCB  specimens  by  Kobayashi,  et  al.  and  Irwin, 
et  al.  were  analyzed  by  this  hybrid  numerical  and  experimental  technique.  Despite 
minor  discrepancies,  the  computed  dynamic  fracture  toughness  and  crack 
arrest  stress  intensity  factors  were  in  reasonable  agreement  with  those  determined 
experimentally.  This  comparative  study  between  different  experimental  setups 
also  indicates  that  the  apparent  differences  in  fracture  dynamic  responses  could 
be  attributed  mainly  to  the  differences  in  material  properties,  bluntness  of  the 
initial  crack  and  specimen  sizes  and  not  to  the  differences  in  experimental  tech- 
niques used. 

INTRODUCTION 

Over  the  past  several  years,  Hahn  et  al.  (1,2,3]  have  been  developing  wedge- 
loaded  single/duplex  double  cantilever  beam  (DCB)  specimens  for  determining  the 
relation  between  dynamic  fracture  toughness,  Kjd,  and  crack  velocity  and  for  mea- 
suring a crack  arrest  stress  intensity  factor,  Kja.  This  specimen  development  was 
accompanied  by  Kannlnen  et  al.'s  comprehensive  one  and  two-dimensional  dynamic 
elastic  analyses  of  the  wedge-loaded  DCB  specimen  (4,5]  with  fixed  grip  loading 
condition.  Later  analytical  developments  by  Kannlnen,  et  al.  Included  the  addition 
of  a test  machine  compliance  in  the  loading  train  for  studying  the  effects  of  ma- 
chine compliance  on  the  dynamic  response  of  a fracturing  DCB  specimen  (6].  The 
dynamic  responses  of  wedge-loaded  DCB  specimens  have  also  been  studied  experimen- 
tally by  dynamic  photoelasticity  (7,8]  and  the  method  of  dynamic  caustics  (9]. 

It  is  not  surprising  that  the  three  series  of  experiments  resulted  in  somewhat 
different  conclusions  regarding  the  dynamic  responses  of  these  DCB  specimens.  The 
results  of  Reference  (7],  for  example,  casts  doubts  on  the  existence  of  a unique 
relation  between  dynamic  fracture  toughness  and  crack  velocity  and  hence  of  a 
crack  arrest  stress  intensity  factor  in  the  Homalite-100  plates  used  for  fracture 
testing.  On  the  other  hand,  a unique  relation  between  dynamic  fracture  toughness 
and  crack  velocity  is  shown  in  Reference  (8]  for  the  same  Homalite-100  material  of 
larger  thickness.  The  crack  arrest  stress  intensity  factor,  Kja,  was  also  found 
to  be  95  percent  of  the  static  fracture  toughness,  Kj  . Post  arrest  stress  inten- 
sity factor  was  also  observed  to  be  slightly  lower  than  Kj  in  agreement  with  the 
concept  of  Kja  based  on  a static  analysis  sometime  after  crack  arrest  (10].  Re- 
cent fracture  testings  of  Aradite-B  specimens  tend  to  confirm  the  above  results 
where  the  crack  arrest  stress  intensity  factor,  Kia,  was  found  to  be  about  equal 
to  the  fracture  toughness.  In  these  experiments,  the  dynamic  stress  intensity 
factors  after  crack  arrest  oscillated  about  the  corresponding  static  value  which 
varied  with  the  crack  velocity  history  (9]  in  apparent  disagreement  with  findings 
of  Reference  (9J. 
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Inherent  In  the  above  widely  varing  conclusions  of  each  series  of  experiments 
was  the  supposition  that  each  result  would  be  generally  applicable  to  any  other  two 
dimensional  dynamic  fracture  problems  regardless  of  sizes,  compliances  of  the  load- 
ing systems  and  static  and  dynamic  material  properties  thus  each  precluding  the 
existence  of  the  other  two  seemingly  contradictory  conclusions.  Before  assessing 
the  possible  variability  in  dynamic  responses  due  to  these  test  parameters,  a 
standard  DCB  specimen  of  common  geometry  and  loading  system  would  have  to  be  ana- 
lyzed by  the  three  groups  of  experimentalists  in  order  to  first  assess  the  experi- 
mental accuracies  of  the  techniques  used.  An  alternate  procedure  would  be  to 
analyze  the  three  different  wedge-loaded  DCB  specimens  with  a common  and  reliable 
analytical  technique.’  The  agreement  or  disagreement  between  the  analytical  and 
experimental  results  could  then  provide  some  insight  into  the  effects  of  specimen 
size  and  material  properties  on  the  dynamic  responses  of  three  different  DCB  spe- 
cimens considered  in  References  [7,8  and  9]. 

The  objective  of  this  paper  is  to  use  such  analytical  procedure  for  a compa- 
rative study  of  the  dynamic  responses  of  one  typical  fracture  test  results  in  each 
of  References  [7,8  and  9]  for  the  purpose  of  deducing  the  effects  of  specimen  geo- 
metries and  material  properties  in  these  three  separate  test  procedures. 

DYNAMIC  FINITE  ELEMENT  ANALYSIS 


The  procedure  used  is  a two-dimensional,  dynamic  finite  element  code,  HONDO 
[11],  which  was  updated  and  modified  for  fracture  dynamic  analysis.*  The  basic 
modifications  consisted  of  algorithms  for  startup  and  for  computing  dynamic  stress 
intensity  factor,  dynamic  energy  release  rate,  fracture  energy,  kinetic  energy  and 
strain  energy  at  each  increment  of  crack  advance. 

In  the  startup  procedure,  the  initial  static  stress  distribution  in  a preload- 
ed structure  prior  to  dynamic  crack  propagation  is  computed.  This  initial  stress 
distribution  must  be  in  complete  static  equilibrium  prior  to  the  initiation  of  a 
dynamic  event.  The  finite  element  breakdown  and  hence  the  initial  stiffness  ma- 
trix used  in  this  preliminary  static  analysis  should  be  identical  to  those  at  the 
initiation  or  at  the  Instant  of  time  t • 0+  in  the  dynamic  analysis.  Close  atten- 
tion must  be  given  to  comoutatlonal  details,  such  as  matchlne  the  2x2  Gaussian 
integration  Dolnts  in  the  orelimlnarv  static  and  subseouent  dvnamic  analvses  in 
order  to  avoid  anv  small  differences  between  the  finite  element  aleorithms  which 
will  be  sensed  as  unbalanced  residual  stresses  and  thus  set  off  oarasltlc  stress 
wave  orooasatlon  in  the  HONDO  II  analysis. 


In  our  oast  dvnamic  finite  element  analvses  of  fracturlne  Homallte-100  dates, 
considerable  oscillations  in  the  calculated  dvnamic  enerev  release  rates  and  hence 
in  the  dvnamic  stress  intensity  factors  were  noted  (12.131.  Althoueh  the  lack  of 
such  oscillations  in  the  corresnondlne  dvnamic  ohotoelastlcitv  results  are  in  oart 
attributable  to  viscous  damoine  in  ohotoelastlc  Dolvmers.  much  of  the  oscillations 
were  thought  to  be  generated  through  the  instant  release  of  crack-tip,  finite  ele- 
ment nodes  during  the  process  of  discrete  crack-tip  advances.  In  order  to  reduce 
the  impulse  stress  waves  generated  by  such  Instantaneous  release  of  a crack-tip 
node,  the  nodal  force  was  reduced  in  equal  increments  which  were  determined  by 
dividing  the  inter-nodal  crack-tip  transit  time  with  the  built-in  finite  time- 
increment  in  HONDO  II.  This  procedure  physically  models  a more  gradual  transit  of 
the  crack-tip  between  two  adjacent  finite  element  nodes.  This  nodal  force  release 
mechanism  is  similar  to  that  developed  by  Keegstra  [14-17]  with  the  exception  that 
the  restraining  nodal  force  is  completely  eliminated  when  the  crack-tip  reaches  the 
adjacent  node.  The  dissipated  energy  during  such  crack  extension  will  be  governed 
by  the  variations  in  nodal  forces  versus  nodal  displacement  relation  during  crack 
extension.  In  general  this  nodal  force  versus  nodal  displacement  relation  is 


* The  updated  finite  element  code  it  referred  to  as  HONDO  II. 
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non-1 inoar  and  will  bo  governed  by  the  dynamic  state  surrounding  the  propagating 
crack  tip  thus  requiring  monitoring  of  nodal  displacement  at  each  incremental 
time  If  the  dissipated  energy  is  used  for  calculating  dynamic  energy  release 
rates.  The  dynamic  stress  Intensity  factor  can  then  be  computed  from  the  dynamic 
energy  release  rate  using  Freund's  relation  [17].  The  generality  of  this  relation 
In  the  presence  of  reflected  stress  waves  In  finite  geometry  was  shown  by  Nilsson 

[18) .  Alternatively,  the  near  field  dynamic  stress  field  as  derived  by  King  et  al. 

[19]  can  be  used  to  calculate  the  dynamic  stress  intensity  factor  directly  from 
the  numerically  obtained  stresses  either  at  the  closest  Gaussian  Integration  point 
or  at  the  center  of  a finite  element  which  shares  the  crack  tip  node. 

The  appropriateness  of  the  above  procedures  for  computing  a dynamic  stress 
Intensity  factor  was  checked  by  analyzing  the  Broberg  problem  [20],  Figure  1 shows 
the  coarse  finite  element  breakdown  used  In  analyzing  a crack  propagating  at  a 
high  speed  of  C/Cj  ■ 0.33  where  C and  Cj  are  the  crack  velocity  and  dilatatlonal 
wave  velocity  In  a steel,  respectively.  The  large  square  finite  element  of  150 
mm  x 150  mm  as  well  as  the  relatively  high  crack  velocity  used  in  this  study  simu- 
lated the  extreme  conditions  experienced  In  another  paper  presented  at  this  Sympo- 
sium and  thus  served  as  an  estimate  of  numerical  errors  Involved  In  the  latter 
[21]. 


Figure  2 shows  the  theoretical  and  computed  crack  opening  displacements  (COD) 
as  the  central  crack  starts  to  extend  from  zero  crack  length  at  constant  rate. 
Despite  the  coarseness  of  the  mesh,  remarkable  agreement  between  the  computed  and 
analytical  CO Ds  at  even  the  first  few  increments  of  crack  extension  Is  noted. 

The  coarseness  of  the  finite  element  mesh  at  the  Initial  phase  of  crack  extension 
suggests  that  the  near  field  COD  equations  from  Reference  [19]  cannot  be  used 
effectively  for  computing  the  dynamic  stress  Intensity  factor,  Kjdyn.  Since  the 
adjacent  Gaussian  Integration  points  and  the  center  of  the  element  was  closer  to 
the  crack  tip,  an  attempt  was  made  to  compute  the  dynamic  stress  Intensity  fac- 
tor, Kjdyn  by  using  the  near  field,  dynamic  state  of  stresses  as  described  in 
Reference  [19].  The  dynamic  stress  Intensity  factors  computed  from  the  normal 
and  devlatorlc  stresses  at  the  nearest  Gaussian  Integration  point,  however,  varied 
as  much  as  AO  percent  from  the  theoretical  values  and  thus  this  procedure  was 
abandoned.  The  dynamic  stress  Intensity  factor  computed  from  the  normal  stress, 
0y„,  at  the  center  of  the  element  as  defined  in  Figure  3 were  more  stable  and 
thus  this  Kjdyn  was  compared  against  the  theoretical  solution  as  shown  in  Figure 
3.  Note  that  much  of  the  spurious  oscillations  in  the  calculated  dynamic  stress 
Intensity  factors  observed  In  previous  analyses  [12,13]  were  eliminated  by  the 
linearly  Increasing  release  of  nodal  force  while  the  crack  tip  advanced  from  one 
finite  element  node  to  another.  The  Initial  large  overestlmatlon  of  Kjdyn,  as 
shown  in  Figure  3,  could  be  attributed  to  the  Inappropriateness  In  using  a one- 
term  representation  of  the  near  field  dynamic  state  of  stress  when  the  crack  ex- 
tended from  zero  crack  length  to  3 to  4 finite  element  lengths.  However,  re- 
markable agreements  between  computed  and  theoretical  Kjdyn  are  noted  for  longer 
crack  length  where  the  one-term  representation  of  the  near  field  dynamic  state  of 
stress  becomes  Increasingly  valid. 

Although  the  above  results  Indicate  the  need  for  finer  element  breakdown  at 
the  Initial  phase  of  the  Broberg  problem,  such  fine  element  breakdown  for  calcu- 
lating Kjdyn  from  the  mid-element  stress  may  not  be  always  practical,  since  the 
time  increment  In  dynamic  finite  element  analysis  Is  governed  by  the  size  of 
Its  smallest  element.  The  strain  energy  release  rate  procedure  of  calculating 
static  stress  Intensity  factors  from  the  results  of  finite  element  analysis,  on 
the  other  hand,  consistently  provided  accurate  static  stress  Intensity  factors 
with  relatively  coarse  meshes  and  thus  the  related  dynamic  energy  release  rate 
procedure  was  used  to  compute  Kjdyn  for  the  same  Broberg  problem.  As  shown  in 
Figure  3,  notable  Improvement  in  the  accuracy  of  K,dyn  at  the  time  of  the  first 
Increment  of  crack  propagation  was  made  but  the  Kjdyn  after  3 to  A Incremental 
crack  extensions  was  not  as  accurate  as  the  Kjdyn  computed  directly  from  the 
mid-element  stress.  Nevertheless,  the  proven  accuracy  of  the  energy  release 
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rate  procedure  in  static  analy- 
sis and  its  reasonable  accuracy 
in  computing  Kjdyn  with  such 
coarse  mesh  of  Figure  1,  i.e. 

150  mm  square,  at  a high  crack 
velocity  of  C/Cj  * 0.33  lead  us 
to  choose  the  procedure  of  dy- 
namic energy  release  rate  for 
computing  Krdyn  in  our  dynamic 
finite  element  analyses  of  the 
wedge-loaded  DCB  specimens  as 
well  as  the  crack  arrest  test 
specimens  [21], 

WEDGE-LOADED  DCB  SPECIMENS 

The  three  wedge-loaded  DCB 
specimens  which  were  analyzed  by 
the  dynamic  finite  element  code 
described  above  are  shown  in 
Figure  4.  For  convenience  in 
identification,  the  three  spe- 
cimens are  designated  as  KML, 
IDKFE  and  KBW  specimens,  re- 
spectively. The  static  and  dy- 
namic material  properties  as 
determined  by  the  three  groups 
of  investigators  [7,8,9)  are 
shown  in  Table  1.  Although  the 
static  material  properties  were 


Figure  3.  Dynamic  Stress  Intensity  Factors  of  a 
Central  Crack  Expanding  at  a Constant 
Rate  in  a Uniaxial  Stress  Field  of  68.9 
MPa  (10,000  psi). 
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SPECIMEN  K0W  (•] 

Three  Hedge-Loaded  DCB  Specimens 


Tabic  1 - Elastic  Properties  of  Wedge-Loaded  DCB  Specimens 


Spec  linen 
Ident 1 f lea. 

Material 

Static 

Dynamic 

Modulus  of 
Elasticity 

GPa 

Poisson's 

Ratio 

Modulus  of 
Elasticity 
CPa 

Poisson's 

Ratio 

KML  [ 7 ] 

Homal lte-100 

3.72 

0.345 

4.65 

0.345 

IDKEF[8] 

Homalite-100 

3.89 

0.31 

4.82 

0.31 

KBW [ 9 ] 

Araldlte-B 

3.38 

0.33 

3.66 

i 

0.39 

comparable,  the  Araldlte-B  epoxy  showed  lesser  strain  sensitivity  and  higher  sta- 
tic fracture  toughness  than  the  two  Homalite-100  plates.  The  30  to  40  percent 
differences  In  static  and  dynamic  elastic  moduli  In  the  Homalite-100  plates  forc- 
ed the  calculation  to  be  conducted  following  the  procedure  [6]  developed  at  Bat- 
telle's  Columbus  Laboratories.  Basically,  the  procedure  is  to  execute  all  static 
and  dynamic  analyses  by  using  the  static  elastic  modulus  and  then  use  the  dynamic 
static  modulus  when  computing  the  dynamic  stress  intensity  factor  from  the  dynamic 
energy  release  rate.*  Identical  fine  meshes  In  the  three  finite  element  break- 
downs, as  shown  in  Figure  5,  were  used  in  analyzing  all  three  specimens  In  order 
to  minimize  the  numerical  errors  due  to  different  fineness  in  finite  element 
breakdown.  The  crack  positions  versus  time  relations  for  the  three  specimens, 
as  shown  in  Plgure  6,  were  then  used  to  drive  the  crack  at  prescribed  rates 
and  the  dynamic  energy  release  rate , ^ i**)'1’,  and  dynamic  stress  intensity  factors, 
Kj^y®,  were  computed  following  the  procedure  described  above.  It  is  interesting 
to. note  that  the  crack  propagated  comparable  distances  in  all  three  specimens 
and  that  the  crack  velocity  in  the  KML  specimen  was  significantly  higher  than 
those  in  the  IDKFE  and  KBW  specimens. 

RESULTS 

KML  Specimen  [ 7 J 

A state  of  plane  stress  was  assumed  in  the  numerical  analysis  of  this  rela- 
tively thin  Homalite-100  plate.  The  calculated  and  measured  dynamic  stress  in- 
tensity factors  as  well  as  the  calculated  static  stress  intensity  factor  versus 
crack  position  are  shown  in  Figure  7.  Since  the  loading  pin  displacement  at 
the  onset  of  crack  propagation  was  not  measured  in  this  series  of  experiments, 
the  stress  intensity  factor  for  crack  initiation,  Kq*t  was  estimated  on  the  basis 
of  matching  the  total  dynamic  energy  released  with  the  calculated  total  static 
strain  energy  released  in  this  specimen.  The  resultant  Kq  would  thus  be  under- 
estimated since  no  estimate  of  the  extraneous  dissipated  energy  in  the  specimen 
is  Included  in  this  calculation.  Reasonable  agreement  existed  between  the  com- 
puted and  measured  Kq  throughout  the  crack  propagation  except  for  the  initial 
phase  of  crack  propagation  and  in  the  region  of  momentary  crack  arrest.  The 
isolated  experimental  point  in  the  former  was  ignored  in  this  comparison  due  to 


* The  superposition  procedure  developed  in  the  original  dynamic  finite  element 
analysis  [12,13]  handles  this  strain  sensitivity  problem  by  using  static  elas- 
tic modulus  In  the  static  calculation  and  dynamic  elastic  modulus  in  the  dyna- 
mic analysis. 

**  Note  that  the  subscript  of  I Is  dropped  for  all  plane  stresj  results. 
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Figure  5.  Finite  Element  Breakdowns  of  Three  DCB  Specimens 

the  blurriness  in  the  dynamic  isochromatic  fringes  and  crack  tip  position  which 
could  have  Introduced  large  errors  in  Kp  determination.  The  minor  discrepancies 
between  the  experimental  and  calculated  Kp  in  the  region  of  crack  arrest  can  be 
attributed  to  the  dynamic  finite  element  analysis  which  is  sensitive  to  the  varia- 
tions in  crack  velocities.  Crack  velocities  measurements  in  this  region  were  not 
accurate  due  to  the  discrete  recording  of  the  crack  which  apparently  arrested 
momentarily  before  starting  up  again. 

IDKFE  Specimen  (8) 

Figure  8 shows  the  variations  in  the  calculated  and  measured  dynamic  stress 
intensity  factors  as  well  as  the  calculated  static  stress  intensity  factors. 

Note  that  the  state  of  plane  strain  was  assumed  in  the  static  and  dynamic  ana- 
lyses of  this  specimen,  not  because  this  Homallte-100  specimen  was  thicker  (13 
mm  versus  10  mm),  but  because  the  plane  stress  results  yielded  a lower  Kp  and 
Increased  the  already  existing  discrepancies  between  measured  and  calculated  re- 
sults. Further  study  of  the  data  in  Table  2.14  and  Figure  2.9  in  Reference  [8] 
indicated  that  perhaps  the  recorded  wedge-pin-openlng  displacement  in  this  ex- 
periment could  be  low  thus  providing  a low  Kjq  on  which  the  entire  static  and 
dynamic  calculations  were  based.  If  Kjq  was  underestimated  by  say  twenty  per- 
cent, then  the  calculated  static  and  dynamic  stress  intensity  factor  curves  will 
shift  upward  and  almost  match  the  experimental  dynamic  stress  Intensity  factors. 


KWB  Specimen  [9] 


Figure  7.  Stress  Intensity  Factors  In  s Wedge- 
Loaded  DCB  Specimen  [7] 


Figure  9 shows  the  varia- 
tions in  the  calculated  and 
measured  dynamic  stress  inten- 
sity factors  as  well  as  the 
calculated  static  stress  inten- 
sity factors.  Reasonable  agree- 
ment in  the  calculated  and  mea- 
sured Kp  are  seen,  with  minor 
differences  in  calculated  and 
measured  values  in  the  region 
of  crack  arrest. 

Figure  10  shows  the  calcu- 
lated variations  in  energies 
with  crack  extension.  These 
energy  variations  follow  the 
characteristic  rapid  decrease 
In  strain  energy,  an  increase 
followed  by  a drop  in  kinetic 
energy  and  gradual  Increase  in 
dissipated  fracture  energy 
[1-6],  The  actual  values  differ 
with  those  in  Figure  2.11  of 
Reference  [6],  particularly  in 
the  former  two  energies.  Part 
of  these  discrepancies  could  be 
attributed  to  the  one-dlmen- 
aional  analysis  used  in  the 
Bat  telle  code  which  would  un- 
derestimate the  strain  energy 
and  hence  the  fracture  energy 
computed  from  energy  balance. 


DISCUSSIONS 


Figure  8. 


Stress  Intensity  Factors  in  a Wedge- 
Loaded  DCB  Specimen  [8]. 


Calculated  and  measured  dy- 
namic stress  intensity  factors  in 
KML  and  KBW  wedge-loaded  DCB  spe- 
cimens agreed  reasonably  well  and 
there  is  reason  to  speculate  that 
similar  agreement  would  have  been 
obtained  in  the  IDKFE  specimen. 

The  dynamic  finite  element  ana- 
lysis reproduced  the  oscillations 
in  Kp  in  the  KML  specimen  as  well 
as  the  relatively  uniform  Kjp  in 
the  IDKFF.  and  KD  in  the  KBW  spe- 
cimen. The  oscillations  in  Kp 
in  the  KML  specimen  could  be  at- 
tributed to  the  smallness  in 
specimen  size,  as  shown  in  Figure 
4,  which  would  generate  higher 
interaction  between  the  reflected 
stress  waves  and  the  propagating 
crack  tip.  This  large  stress 
wave  effect  was  further  augmented 
by  the  high  Kq  value  necessary  to 
drive  the  crack  approximately 
the  same  distance  as  in  the  other 
two  IDKFE  and  KBW  specimens. 

The  computed  overshoot  in  Kp  im- 
mediately after  crack  propagation 
could  also  be  attributed  to  the 
large  stress  wave  effect  in  the 
KML  specimen.  The  lower  crack 
initiation  stress  intensity  fac- 
tors in  the  IDKFE  and  KBW  speci- 
mens combined  with  the  much  longer 


specimen  sizes  obviously  diminished  the  stress  wave  effect  as  shown  by  the  lack  of 
oscillations  in  the  experimental  and  numerically  determined  dynamic  stress  inten- 
sity factors. 

The  gradual  deceleration  crack  speed  prior  to  crack  arrest  and  thus  the  exis- 
tence of  a distinct  crack  arrest  stress  intensity  factor,  KIa,  are  noted  in  the 
IDKFE  and  KBW  specimens.  The  high  static  fracture  toughness,  Kjc,  of  Araldite-B 
could  be  responsible  for  the  closeness  in  Kja  and  Kjc  in  the  KBW  specimen  as  the 
crack  slows  down  to  an  arrest.  Kja  in  the  KML  specimen  is  less  distinct,  possibly 
due  to  the  lack  of  experimental  data  at  finer  time  increments  during  the  period 
of  momentary  crack  arrest.  Again  the  difference  between  the  crack  arrest  charac- 
teristics could  be  attributed  to  the  differences  in  ((^specimen  sizes  and  the 
associated  stress  wave  effects. 

The  calculated  and  measured  dynamic  arrest  stress  intensity  factors  of  the 
three  specimens  were  always  lower  than  the  corresponding  measured  fracture  tough- 
nesses, Kj  , and  higher  than  the  corresponding  static  stress  intensity  factor.  The 
variability  in  the  latter  static  stress  intensity  factor,  as  noted  in  Figures  7, 

8 and  9,  probably  exclude  this  value  as  material  property  related  to  crack  arrest. 

CONCLUSIONS 


The  updated  HONDO  II  dynamic  finite  element  code  with  incremental  release 
of  crack  tip  nodal  force  has  been  shown  to  be  a reliable  procedure  in  analyzing 
fracture  dynamic  problems. 


I 


Figure  9.  Stress  Intensity 
Factors  In  a Wedge-Loaded  DCB 
Specimen  [9]. 


Figure  10.  Energies  In  a 
Wedge-Loaded  DCB  Specimen  (9). 


This  code  succersfully  duplicated  the  experimental ly  determined  dynamic 
fracture  toughness  in  two  of  the  three  fracturing  wedge-loaded  DCB  specimens  and 
showed  that  the  apparent  differences  in  fracture  dynamic  responses  could  be  at- 
tributed mainly  to  the  differences  in  material  properties,  bluntness  of  the  ini- 
tial crack  and  specimen  sizes  and  not  to  the  differences  in  experimental  tech- 
niques used. 

ACKNOWLEDGEMENT 

The  results  of  this  investigation  were  obtained  in  a research  contract 
funded  by  the  Office  of  Naval  Research  under  Contract  No.  N00014-76-C-0060, 

NR  064-478.  The  authors  wish  to  acknowledge  the  support  and  encouragement  of 
Drs.  N.R.  Perrone  and  D.  Mulville  of  ONR  during  the  course  of  this  investigation. 
The  authors  also  wish  to  acknowledge  the  discussions  with  Professor  W.L.  Fourney, 
University  of  Maryland  and  with  Dr.  J.F.  Kalthoff,  Institut  fur  Festk&'rpermech- 
anik. 

REFERENCES 

1.  HAHN,  G.T.,  HOAGLAND,  R.G.,  KANNINEN,  M.F.  and  ROSENFIELD,  A.R.  - Pilot 
Study  of  Fracture  Arrest  Capabilities  of  A533B  Steel.  Cracks  and  Fractuie, 
ASTM  STP  601,  pp.  209-233,  June  1976. 

2.  HOAGLAND,  R.G.,  CEHLEN,  P.C.,  ROSENFIELD,  A.F.  and  HAHN,  G.T.  - Characteris- 
tics of  a Run-Arrest  Segment  of  Crack  Extension.  Fast  Fracture  and  Crack 
Arrest.  ASTM  STP  627,  pp.  203-207,  July  1977. 

3.  HAHN,  G.T.,  HOAGLAND,  R.G . and  ROSENFIELD,  A.R.  - A Fracture  Mechanics 
Practice  for  Crack  Arrest.  Trans,  of  the  4th  Int'l,  Conf.  on  Structural 
Mechanics  in  Reactor  Technology,  CECA,  CEE,  CEEA  Luxembourg,  Paper  G 1/6, 

1977. 

4.  KANNINEN,  M.F.  - A Dynamic  Analysis  of  Unstable  Crack  Propagation  and  Arrest 
in  the  DCB  Test  Specimen.  Int'l.  J.  of  Fracture,  Vol.  10,  No.  3,  pp.  415- 
431,  September  1974. 

5.  KANNINEN,  M.F.,  POPELAR,  C.  and  GEHLEN,  R.C.  - Dynamic  Analysis  of  Crack 
Propagation  in  the  DCB  Specimen.  Fast  Fracture  and  Crack  Arrest,  ASTM  STP 
627,  pp.  19-38,  July  1977. 

6.  HAHN,  G.T. , GEHLEN,  R.C.,  HOAGLAND,  R.G.,  MARSHALL,  C.W.,  KANNINEN,  M.F., 
POPELAR,  C.  and  ROSENFIELD,  A.F.  - Critical  Experiments,  Measurements  and 
Analyses  to  Establish  a Crack  Arrest  Methodology  for  Nuclear  Pressure  Ves- 
sel Steels.  Task  62,  Second  Annual  Report,  Battelle  Columbus  Laboratories 
BMI-1959,  October  1976. 

7.  KOBAYASHI,  A.S.,  MALL,  S.  and  LEE,  M.H.  - Fracture  Dynamics  of  Wedge-Loaded 
DCB  Specimen.  Cracks  and  Fracture,  ASTM  STP  601,  pp.  274-290,  June  1976. 

8.  IRWIN,  G.R.,  DALLY,  J.W.,  KOBAYASHI,  T. , FOURNEY,  W.L.  and  ETHERIDGE,  J.M. 

- A Photoelastic  Characterization  of  Dynamic  Fracture.  University  of  Maryland 
Report  prepared  for  the  U.S.  Nuclear  Regulatory  Commission,  NUREC-0072, 

NRC-5,  December  1976. 

9.  KALTHOFF,  J.,  BEINERT,  J.  and  WINKLER,  S.  - Measurements  of  Dynamic  Stress 
Intensity  Factors  for  Fast  Running  and  Arresting  Cracks  in  Double-Cantilever- 
Beam  Specimens.  Fast  Fracture  and  Crack  Arrest,  ASTM  STP  627,  pp.  161-176, 
July,  1977. 

10.  CROSLEY,  R.P.  and  RIPLINC,  E.J.  - Characteristics  of  a Run-Arrest  Segment  of 
Crack  Extension.  Fast  Fracture  and  Crack  Arrest.  ASTM  STP  627,  pp.  203-227, 
July  1977. 

11.  KEY,  S.W.  - HONDO,  A Finite  Element  Computer  Program  for  the  Large  Defor- 
mation Dynamic  Responses  of  Axisymmetric  Solids,  Sandla  Laboratories. 

12.  KOBAYASHI,  A.S.,  EMERY,  A.F.  and  MALL,  S.  - Dynamic  Finite  Element  and 
Dynamic  Photoelastic  Analyses  of  Crack  Arrest  in  Homalite-100  Plates,  Fast 
Fracture  and  Crack  Arrest.  ASTM  STP  627,  pp.  95-108,  July  1977. 

13.  KOBAYASHI,  A.S.,  EMERY,  A.F.  and  MALL,  S.  - Dynamic  Finite  Element  and 


Dynamic  Photoelastic  Analysis  of  Two  Fracturing  Homalite-100  Plates.  Experi- 
mental Mechanics.  Vol.  16,  No.  9,  pp.  321-328,  September  1976. 

14.  KEEGSTRA,  P.N.R.  - A Transient  Finite  Element  Crack  Propagation  Model  for 
Nuclear  Reactor  Pressure  Vessel  Steels.  J.  Inst.  Nucl.  Engrs..  Vol.  17.  No. 

4,  pp.  89-96,  1976. 

15.  KEEGSTRA,  P.N.R. , HEAD,  J.L.  and  TURNER,  C.E.  - A Transient  Finite  Element 
Analysis  of  Unstable  Crack  Propagation  in  Some  2-Dlmenslonal  Geometries. 

Proc.  of  the  4th  Int'l.  Conf.  on  Fracture.  University  of  Waterloo  Press, 

Vol.  3,  pp.  515-522,  1977. 

16.  KEEGSTRA,  P.N.R.,  HEAD,  J.L.  and  TURNER,  C.E.  - The  Interpretation  of  the 
Instrumented  Charpy  Test.  Trans,  of  the  4th  Int'l.  Conf.  on  Structural 
Mechanics  in  Reactor  Technology,  Vol.  G.,  CECA,  CEE,  CEEA  Luxembourg,  paper 
cf  4/7,  1977. 

17.  FREUND,  L.B.  - Crack  Propagation  in  an  Elastic  Solid  Subjected  to  Ceneral 
Loading  - II  non-Uniform  Rate  of  Extension.  Journal  of  Mechanics  and  Physics 
of  Solids,  Vol.  20,  1972,  pp.  141-152. 

18.  NILSSON,  F.  - A Note  on  the  Stress  Singularity  at  a Non-Uniformly  Moving 
Crack  Tip.  Journal  of  Elasticity.  Vol.  4,  No.  1,  March  1974,  pp.  73-75. 

19.  KING,  W.W.,  MALLUCK,  J.F.,  ABERSON,  J.A.  and  ANDERSON,  J.M.  - Application  of 
Running  Crack  Eigenfunction  to  Finite  Element  Simulation  of  Crack  Propagation. 
Mechanics  Research  Communication,  Vol.  3,  No.  3,  pp.  197-202,  1976. 

20.  BORBERG,  K.B.  - The  Propagation  of  a Brittle  Crack.  Arkiv  fur  Fysik,  Vol.  18, 
pp.  159-198,  1960. 

21.  KANAZAWA,  T.,  KOBAYASHI,  A.S.,  MACHIDA,  S.  and  URABE,  Y.  - Fracture  Dynamic 
Analysis  of  Crack  Arrest  Test  Specimens.  To  be  published  in  the  Proc.  of 
this  Symposium. 


PAD’  ' i W'tr(RN*WNT 

‘•'•bi.tr.t.,,  m i.jtj.n  Actl.HlM 

0»fi  r «f  Naval  «*SMr<h 

••  tm*-i  f !►»  \j v > 

*r'  V /?;-?; 

•'ttn  iati  4;4 

an 

POO 


Olr*K  tor 
• . . , f y 

1 1 .i  rb-vn.ci  , h 

1 - urr..  ■ 

|r«*et 

022 IC 

Diret  tor 

f,<f ' ■ '<***  ' °.C<.rar  h 

N.H  Snytr  i i 

srk  Street 

• • • 

fObOS 

Oirec tor 

0M»  ••  • N* 

*al  Research 

’“■**  ’ -r*  <•«  ‘iff.  ,, 

1 Rmedway 

1 1 

'ie»  tnrfc, 

ID'VM 

• ” Branch  nfflce 

' * • ‘ Sf  e»r  Strot 
Pa-.«(iena,  (.A  9110ft 

D Ire.-  tor 

**  * Naval  Research 

,'*"CU<6  Area  'tfficc 

’nr  . ;,r  P1a/a,  Su,tf  40) 

San  f r«n<  isco,  CA  Mio ? 

*r'earch  tal-crator*  *61 
tod*  ?6?;  ' * # 

Washington.  D C.  20375 

••’f  -.p  Doowentatlom  r»n>»  /ij) 
>t4t»on  ' 

A'r.arxJri*.  vn  ??J14 
NAVY 

V'.a)  ^-sear-  h laturat^r, 

*'■»  hlngton,  0 c . J0JJ5 

»ttn  Code  MOO 
8410 
M30 
8440 
MO0 
Cl  90 
6 MO 

" m-  Ot  rarCn  t Hey  Center 
' U **«>'•'  Shipyard 
Portsmouth,  V«  ?37i9 
Attn  Op  t Palmer.  Code  177 


havy  (Continue) 

“T*  i.  '"'P  *««rc» 

and  Dtvylnpufnt  Certer 
Annapolis.  HO  2J40. 

Attn  coda  /74c 
?R 
?8I 

J!;S  ^t*1  Capons  Center 
li’na  lake,  CA  9155s. 

Attn  Code  406/ 

4520 

C.omuniinq  Office 

C*.  ul*’  C,,"  laboratory 

9ort  M0eneme,  rj  9J04J 

9ay.u  Surface  Weapons  Center 
9h>  *e  Oak 

Stl*e.  Spring,  *n  /09)0 
Attn  code  WP-lo 
WA-’o 

Tf;  hni<  si  Otrer  tor 
pava)  0 <ean  Systems  Center 
■*n  Diego.  CA  9?»52 

Supervisor  of  Shipbuilding 
U.S  Hj*y 

Newport  News.  VA  ?J607 

wiiiMjsrn&isr1 

P.C  Bo«  8317 
Orlando.  Fl  32806 

Chief  of  Naval  Operations 
department  0f  the  Navy 

Washington , oc  ?0350 
Attn  rode  OP-098 

Strategic  Systems  Project  Office 

Ueuar tment  of  the  Nary 
‘Washington,  DC  20376 
Attn:  NSP-200 
Naval  Air  Systems  Co-wand 
department  of  the  Navy 
Washington,  r>C  20361 

CM'  ;**roioaca  | Structgrt.) 

J'orbolcal  ll»r.ry| 

3?08  (Structures) 

Naval  Air  Development  Center 
Director,  Aerospace  Mechanics 
Warminster,  PA  18974 
U.S.  Nava)  Academy 
Engineering  Department 
Annapolis,  IP  71402 


NASA  ('Continued/ 

SrlenUfir  A Technical  Information  Factll 
NASA  Representative  (S-At/DL! 

P 0.  Rom  5700 
Bethesda.  HP  20014 

Other  Qnvtn— at  a.  ttvittei 

'"onvriandan  f 

Chief.  Testing  A Development  Division 
OS  Coast  Guard 
1300  f Street.  NW 
Washington,  DC  20226 
Technical  Director 
War  me  Corps  Development  and 
education  Command 
Quant  1 co,  /A  2?) 34 

Oi rector 

National  Bureau  of  Standards 
Washington.  DC  20034 
Attn:  Mr  B. I . Wilson,  EM  ?19 
br  M Gaus 

National  Science  Foundation 
Envlronr-ental  Research  Division 
Washington,  DC  20550 
library  of  Congress 
Science  and  Technology  Division 
Washington,  DC  20540 
Director 

Defense  Nuclear  Agency 
Washington,  DC  20305 
Attn:  SPSS 

Director  Oefense  Research  S Fnglneerlno 
Technical  tlbrary 
Room  3C128 
The  Pentagon 
Washington,  DC  20301 
Mr.  Jerome  Persh 
Staff  Specialist  for  Materials 
and  Structures 
ODORSE , The  Pentagon 
Room  3D 1 089 
Washington,  DC  20301 

Chief,  Airframe  4 Equipment  drench 
PS -120 

Office  of  Flight  Standards 
Federal  Aviation  Agency 
Washington.  DC  20553 

Chief,  Research  and  Development 
Maritime  Administration 
Washington,  DC  20235 
Plcatlnny  Arsenal 

Plastics  Technical  Evaluation  Center 
Attn  Technical  Information  Center 
Dover,  NJ  07801 


Other  Government  Activities  (Continued) 

Deputy  Chief,  Office  of  Ship  Constructor 
Maritime  Administration 
Washington.  DC  70235 
Attn  Mr  U ( Russo 
National  Academy  of  Science1, 

National  Resea. -ch  Council 
Ship  Hull  Pe-.earef*  Committee 
2101  Constitution  Avenue 
Washington,  Of  20418 
Attn  Mr  4 9 l yt  ?e 
National  Science  Inundation 
Engineer tng  Mechanics  Section 
Division  of  Engineering 
Washington,  DC  20550 
Coewi *r>der  f ie!d  Co-mend 
Defense  Nuclear  Agency 
Sandta  ease 

Albugueroue.  W<  87115 
Atomic  Energy  rcmwission 
Oiv.  of  P-eac  tor  Dev  t.  Tech 
Germantown,  “ary land  '0?f.  1 


**y/  (Continued) 

Naval  facilities  Engineering  'anMid 
200  Store >7  Street 
Ale»andna.  VA  22332 

Attn  Code  03  (Research  4 Ueveiommnt) 

045 

14114  (Technical  library/ 

Na.al  Sea  Systems  r newer.; 
beiiertment  n*  the  Navy 
Washington.  DC  .'016. 

Attn  Code  V (Pesej..h  4 technology) 

I.’  (Ship  S' len  <ng  flivi-.  ion) 
015  (Mechan'cs  | “eterlals) 
**aval  Ship  Engineering  enter 
Department  0t  'he  «u.r 

Washington.  DC  20367 
Attn  Code  6106G 
6*14 
61200 
6128 
6129 

fommanding  Officer  and  Director 
David  W Taylor  Naval  Ship  Research 
and  Development  Tenter 
Bethesda.  *D  20034 
Attn . Code  042 
17 

172 

173 

174 

l§00 
110?  1 
1900 
IfOl 
INI 
I960 
1962 

Naval  tMdeovater  Systems  Center 
"•wp ort.  »f  02440 
Attn:  Dr  R Tralftor 

Neva)  Surface  Weapont  Center 
DRhlgrtn  Laboratory 
Dahlgren,  VA  22448 
Attn:  Code  06-20 
DC- 30 

Teehnfee)  Director 
•Nre  Island  Nava)  Shipyard 
veilejo,  CA  94592 

Camanding  Officer  (?) 

0-5  Anhy  Research  Office 
70  8oi  12211 

••search  Trianel*  hark,  HC  22709 
Attn-  Mr  J.J  Murray,  CRD-AA-IP 


Army  I Continued) 


Water  , l iet  Ars.-nal 
MAGGS  Research  Center 
Wetervl iet , NT  12)89 
Attn  Olrector  of  Research 

U S Army  Materials  and  Mechanics 
Research  Center 
Watertown,  HA  02172 

Attn-  Dr  R Shea,  DRXMR-T 


B "Mint  Research  A Dev  CenU 

Redstone  Scientific  Information  Center 
Chief.  Document  Section 
Redstone  Arsenal,  At  15809 


Army  Research  and 
Fort  Bel  voir,  VA 


Development  Center 

22060 


•t  r Force 


Commander  MA DO 

Wrtght-Patferion  Air  Force  Base 
Dsyton,  OH  45433 
A'tn  Code  MMMDC 

AF701  (FOOS) 

Strut  turn  Otviiion 
AFIC  (MCEEA) 


.....  . . .-VJ  -»in.nii5  uroup 

»lr  -ort.  iMtltut,  of  TKbnolojy 
Wright 'Patterson  Air  Force  Base 
Olyfon,  OH  454J3 


tnier.  civil  Engineering  Branch 
■IRC.  Research  Division 
• ir  For-ie  Weapons  laboratory 
Mrtland  Air  Force  Base 
Albuguergue.  m 871)7 


Air  Force  Office  of  Scientific  Reseerrh 
Bolling  Air  Force  Bete 
Wtshlngton.  DC  20332 
Attn:  Mechanics  Division 
Department  of  the  Air  Force 
Air  Uni  vanity  library 
Manwell  Air  force  Bata 
Hontgomery.  At  36112 


National  Aeronautics  I Space  Admlnlstratt 
Structures  Research  Olvtslon 
Langley  Res earch  Center 
lRf»9»Ry  Station 
HR»pton,  VA  23365 

Netlonal  Aeronautics  » Spec#  A4»in1itratt 
Associate  Administrator  for  Advanced 
Research  and  Technology 
Washington.  DC  12546 


PART  2 - CONTRACTORS  AND  OTHER  TECHNICAL  COLLABORATORS 


Universities 
Dr. J.  Tinsley  Oden 
University  of  Te*as  at  Austin 
345  Engineering  Science  Building 
Austin,  TX  7871? 

Professor  Julius  Milrlowlte 
California  Institute  of  Technology 
Olv.  of  Engineering  * Applied  Sciences 
Pasadena.  CA  91109 


Or.  Harold  Uebowltx,  Dean 
School  of  Engineering  4 Applied  Science 
George  Washington  University 
Washington,  DC  2005? 

Professor  EH  Sternberg 
California  Institute  of  Technology 
Dlv.  of  Engineering  6 Applied  Science 
Pasadena.  CA  91109 


Professor  Paul  M.  Naghdi 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley-.  CA  94? 20 

Professor  P.S  Symonds 
Brown  University 
Division  of  Engineering 
Providence.  RI  0291? 

Professor  A, J.  Ourellf 
Oakland  University 
School  of  Engineering 
Rochester.  Ml  48063 
Professor  F.i.  OINagglo 
Columbia  University 
Department  of  Civil  Engineering 
Hen  NOrk.  NY  10027 


rruressor  norman  Jones 

Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge.  MA  02 119 


Professor  E.J  Skudriyfc 
Pennsylvania  State  University 
Applied  Research  laboratory 
Department  of  Physics 
State  College,  RA  16801 

Professor  J Kempner 

Polytechnic  Institute  of  New  York 

Dept  of  Aerospece  Ingrg  I Applied  Wech 

)3J  Jay  Street 

Brooklyn,  NY  1 1?Q1 


Professor  j Flosner 

Polytechnic  In-.titute  of  New  York 

Orp»  of  Aerospace  Engrg  s Applied  Me«han.is 

333  Jay  Street 

Brooklyn,  h*  lj?0) 


Professor  I#  Schapery 
Texas  AAM  University 
Department  of  C<v<!  Engineering 
College  Station.  TX  77843 


Professor  Walter  D.  Pilkey 
University  0f  Virginia 
Research  laboratories  for  the  Engineering 
Sciences 


Scbool  of  Fngtn.,r1ng  ,o«  Science 

Charlottesvll le.  VA  2290 1 


Professor  K.Q.  Wlllmert 
Clarkson  College  of  Technology 
Department  of  Mechanical  Engineering 
Potsdam.  NY  13676 


university  0T  Maryland 
Aerospace  Engineering  Department 
College  Park,  MO  2074? 

Or  Walter  f.  Haister 
Tevas  A AM  University 
Aerospace  Engineering  Department 
College  Station.  Tx  77843 


Dr.  BV  Berger 
University  of  *rylang 
Department  of  Mechanical  Engineering 
College  Park.  MO  ?0?4? 


Dr.  I A Schmlt 
Uhlvers tty  of  California 
School  of  Engineering  | Applied  5c 
los  Angeles.  CA  90024 


tence 


university  of  Arlyona 

Dept  of  Aerospace  I Mechanical  Engineering 
Tucson,  A I 85721 


Dr.  S.J.  Fenves 

Carnegie -Me  11  on  University 

Department  of  Civil  Engineering 

Schenley  Park 

Pittsburgh,  RA  15213 

Dr  Ronald  l.  Huston 

Dept  of  Engineering  Analysis 

University  of  Cincinnati 

Cincinnati.  OH  45221 


Universities  (Continued) 

Professor  G.C.M.  Sih 
Lehigh  University 
Institute  of  Fracture  and 
Sol  id  Mechanics 
Bethlehem,  PA  18015 

Professor  Albert  S.  Kohayashi 
University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  WA  981 90 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
State  University 
Dept,  of  Engineering  Mechanics 
Blacksburq.  VA  ?4061 

Professor  A.C.  Eringen 

Dept,  of  Aerospace  A Mech.  Sciences 

Princeton  University 

Princeton,  NO  08040 

Professor  F.H.  Lee 
Stanford  University 
Oiv.  of  Engineering  Mechanics 
Stanford,  CA  94305 

Professor  Albert  I.  King 
Wayne  State  University 
Biomechanics  Research  Center 
Detroit.  MI  48202 

Dr.  V.R.  Hodgson 
Wayne  State  University 
School  of  Medicine 
Detroit.  MI  48202 

Dean  B.A.  Boley 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  IL  60201 

Professor  P.G.  Hodge,  Jr. 

University  of  Minnesota 
Dept,  of  Aerospace  Engineering 
and  Mechanics 
Minneapolis,  MN  55455 

Dr.  DC.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbana.  IL  61801 

Professor  N.M.  Newmark 
University  of  Illinois 
Dept,  of  Civil  Engineering 
Urbana,  IL  61R03 

Professor  E.  Reissner 
University  of  California,  San  Diego 
Dept,  of  Applied  Mechanics 
La  Jolla,  CA  92037 


Universities  (Continued) 

Professor  Kent  R.  Wilson 
University  of  California,  San  Diego 
Department  of  Chemistry 
La  Jolla,  CA  92093 

Professor  B.K.  Donaldson 
University  of  Maryland 
Aerospace  Engineering  Department 
College  Park,  MD  20742 

Professor  Joseph  A.  Clark 
Catholic  University  of  America 
Dept,  of  Mechanical  Engineering 
Washington.  DC  20064 

Professor  T.C.  Huang 
University  of  W1 scons  In-Madison 
Dept,  of  Engineering  Mechanics 
Madison.  WI  53706 

Dr.  Samuel  B.  Batdorf 

University  of  California 

School  of  Engineering  6 Applied  Science 

Los  Angeles,  CA  90024 

Industry  and  Research  Institutes 
U.S.  Naval  Postgraduate  School 
Library 
Code  0384 

Monterey.  CA  93940 

Webb  Institute  of  Naval  Architecture 
Attn:  Librarian 
Crescent  Beach  Road,  Glen  Cove 
Long  Island,  NY  11142 


Professor  William  A.  Nash 

University  of  Massachusetts 

Dept,  of  Mechanics  A Aerospace  Engineering 

Amherst,  MA  01002 

Professor  G.  Herrmann 
Stanford  University 
Department  of  Applied  Mechanics 
Stanford,  CA  94305 

Professor  J.D.  Achenbach 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  IL  60201 

Professor  G.R.  Irwin 
University  of  Maryland 
Dept,  of  Mechanical  Engineering 
College  Park,  MD  20742 

Professor  SB.  Dong 
University  of  California 
Department  of  Mechanics 
Los  Angeles,  CA  90024 

Professor  Burt  Paul 
University  of  Pennsylvania 
Towne  School  of  Civil  and 
Mechanical  Engineering 
Philadelphia,  PA  19104 

°rofessor  H.W.  Liu 
Syracuse  University 

Dept,  of  Chemical  Engineering  A Metallurgy 
Syracuse,  NY  13210 

Professor  S.  Bodner 
Technion  RAD  Foundation 
Haifa,  Israel 

Professor  Werner  Goldsmith 
University  of  California 
Dept,  of  Mechanical  Engineering 
Berkeley,  CA  94720 

Professor  R.S.  Rlvlln 
Lehigh  University 

Center  for  the  Application  of  Mathematics 
Bethlehem,  PA  18015 

Professor  F.A.  Cozzarelli 
State  University  of  New  York  at  Buffalo 
Div.  of  Interdisciplinary  Studies 
Karr  Parker  Engineering  Building 
Chemistry  Road 
Buffalo.  NY  14214 

Professor  Joseph  L.  Rose 
Drexel  University 

Dept,  of  Mechanical  Engineering  A Mechanics 
Philadelphia,  PA  19104 


line  lass i f i ed 


Itrunev  Cl  » /I»|f  ATION  or  Twr,  Pl^f  />•»•  rntorod) 


REPORT  DOCUMENTATION  PAGE 


2 GOVT  ACCESSION  NO 


C>\  A Numerical  Dynamic  Fracture  Analyses  o 
Three  Wedge- Loaded  DCB  Specimens  < 


>-  R^»wOR-») — - > 

\.S.  Aobavashi , S./Mall,  Y./jrabc  and 
A.iF.Mmerv 


• y^:mcr\ 


forming  organization  name  and  aooress 

University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  Washington  98195 

II  CONTROLLING  OFFICE  NAME  AnO  ADDRESS 

Office  of  Naval  Research 
Arlington,  Virginia 

r*  monitoring  AGENCY  name  * ADDRESSfU 


KKAl)  INSTRUCTIONS 
flKKOKi:  COMI'I.I  TINC,  FORM 


I RCCiriEn  t s c at  alog  numbe  h 

TR  No.  29 


». — tmpC  am  ilcnul  a par. on  COVERED 


. j r\ 

Interim /Repygt  / 


$ PCHFOHi/nG  ONG. 


gomthact  %m  £*(•) 

NfV^/)  1 l -7()-c-M(it^y 
NR  064-475' 


10  PROGRAM  ELtMFNT,  PftOJFCT,  TASK 
AREA  a WORK  UNIT  NUMBERS 


t?  RCRtfR*  pure  ■—  ■— 1 
0c  tutor  k$77 


i: 


— - ■ 


is  security  class,  (ot  tnio  topon) 

unclassified 


1S«  OECL  ASSIFIC  ATION/OOWNGRAOiNG 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  fol  this  Report) 


unlimited 


muBUTiON  sTAT:-:.v.':Nf  a 1 


ved  for  public  release; 
Distribution  Unli.T.i  ed 


17  DISTRIBUTION  STATEMENT  (of  iho  oboirmct  ontorod  In  Block  20.  II  dltloronl  Irom  Report) 


18  SUPPL EMENT ARY  NOTES 


19  KEY  WOROS  (Contlnuo  on  rovoroo  aldo  ll  nocooiory  mnd  Idonllty  by  block  numboe) 


Fracture  Mechanics 
Crack  Propagation 
Crack  Arrest 


Impact  — 

Dynamic  Photoelasticity 


ABSTRACT  (Contlnuo  on  rmvoeoo  tldo  II  nocoooory  ond  Idonllly  by  block  numbor) 

A dynamic  finite  element  code  is  used  to  compute  the  dynamic  fracture  toughness 
and  crack  arrest  stress  intensity  factor  from  experimentally  determined  crack 
velocities  in  three  fracturing  wedge-loaded  double  cantilever  beam  (DCB)  speci- 
mens. One  experiment  involving  an  Aradite-B  DCB  specimen  by  Kalthoff,  et  al.,and 
two  experiments  involving  llomalite-100  DCB  specimens  by  Kobayashi,  et  al.  and 
Irwin,  et  al.  were  analyzed  by  this  hybrid  numerical  and  experimental  technique. 
Despite  minor  discrepancies,  the  computed  dynamic  fracture  toughness  and  crack 

I i' mil  liiiiud  mi  ruv.u-L.'l  ~ -.i~ — ■■ — — — 


DD  , 1473 


COITION  OF  I MOV  •*  IS  OSIOLCTC 
S/N  0 10  1- n I 4-  tto  1 I 


J-0 14-  660 1 I , . , / — 

• '/0  3 W 


Unc lass ified 


■ - ■ - — — S*4 

RITV  CL  A»tlf  ICATION  OF  THIS  >‘*GE  Polo  Knlotod) ^ 


w 

> 


ip.  (Continued) 


arrest  stress  intensity  factors  were  in  reasonable  agreement  with  those 
determined  experimentally.  This  comparative  study  between  different 
experimental  setups  also  indicates  that  the  apparent  differences  in 
fracture  dynamic  responses  could  be  attributed  mainly  to  the  differences 
in  material  properties,  bluntness  of  the  initial  crack  and  specimen  sizes 
and  not  to  the  differences  in  experimental  techniques  used. 
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